October, 1998 
hep-ph/9810336 



Possible new physics signals 

in b — > S7 and b — » s/ + /~ 



KAIST-TH 98/15 
KIAS-P98021 
SNUTP 98-078 



OS 
On 
On 

o3 



(N 

(N 
> 

m 
m 
o 

00 
ON. 

Oh: 



X 



Yeong Gyun Kim and P. Ko 
Department of Physics, KAIST 
Taejon 305-701, Korea 

Jae Sik Lee 
School of Physics, KIAS 
Seoul 130-012, Korea 

Abstract 

We consider possible new physics contributions to b — > sl + l~ assuming the 
new physics modifies (chromo) magnetic and electric form factors in b — > S7* 
and b — > sg with the same chirality structure as in the standard model. 
Parametrizing the new physics effects on b — > sj* and b — > sg in terms of four 
real parameters, one finds that there are enough region of parameter space 
in which the measured branching ratio for B — > X s j can be accomodated, 
and the predicted CP violation effect could be as large as ~ 30%. More- 
over, the branching ratio and the forward-backward asymmetry of a lepton 
in B — > X s l + l~ and the tau polarization asymmetry in B — ► X s t + t~ can 
be deviated from the SM predictions by a factor of ~ 2, which can be acces- 
sible at B factories. We also discuss these observables in a specific class of 
supersymmetric models with gluino-mediated flavor changing neutral current 
(FCNC). 
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I. INTRODUCTION 



The missions of B— factories under constructions are (i) to test the CP violation in the 
Standard Model (SM) a la Kobayashi-Maskawa scheme 0, and (ii) to find out any new flavor 
violation and especially new source of CP violation beyond the KM phase in the SM with 
three generations. The latter is well motivated by the fact that the KM phase in the SM may 
not be enough to generate the baryon number asymmetry in the universe. In terms of physics 
view point, the second mission seems more exciting one, since it could uncover a veil beyond 
the SM and provide an ingredient that is necessary to explain baryon number asymmetry 
of the universe. Then, one has to seek for a possible signal of new physics in rare decays of 
B— mesons and CP violation therein. One could choose his/her own favorite models to work 
out the consequences of such model to the physics issues that could be investigated at B 
factories. Or one could work in the effective field theory framework, in a manner as much as 
model-independent as possible. In the following, we choose the second avenue to study the 
possible signals of new physics that could be studied in detail at B factories. Then we give 
explicit examples (that satisfy our assumptions made in the model independent analysis) in 
super symmetric (SUSY) models with gluino-mediated b — > s transition. 

If one considers the SM as an effective field theory (EFT) of more fundamental theories 
below the scale A, the new physics effects will manifest themselves in higher dimensional 
operators (dim [O] > 5) that are invariant under the SM gauge group. Several groups have 
made a list of dimension-5 and dimension-6 operators in the last decade 0. Assuming the 
lepton and baryon number conservations, there are about 80 operators that are independent 
with each other. It would be formidable to consider all of such operators at once, even if 
we are interested in their effects in B physics. However, if we restrict to b — > 57, only two 
operators become relevant : 

v v 
aL ^2 sl<V&a F^, and a R —s R a^b L , (1) 

after the electroweak (EW) symmetry breaking (v is the Higgs vacuum expectation value). 
Here o^r's are dimensionless coefficients. Thus the above operators can be recasted into 
the following form [] : 

H eS (b - 57) = [C 7L 7L + C m O m ] , (2) 

where At = V* s V t b{= —AX 2 in the Wolfenstein parametrization ||]) and 

7L = ^ m b sla^b R F, v . (3) 

The operator 7R is obtained from 7 l by the exchange (L <-> R). Similarly one can expect 
a new physics contribution to b — > sg : 

H cS (b - sg) = [C 8L 8L + C 8R O SR ] , (4) 



L We follow the convention of Ref. H 
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where 



Osl = ^ m b sta^T a a / R G%, (5) 



and 8 r is obtained from 8i by the exchange (L ^ R). These two processes b — > 37 and 
b — > sg are unique in the sense that they are described in terms of only two independent 
operators 07(s)l and 07(8)r whatever new physics there are. This fact makes it easy to study 
these decays in a model indepdent manner ||. 

The SM predictions for the Cr t s at the Mw scale are (in the limit m s = 0) 

CfM (M w ) « -0.22, 
C^(M W ) = 0, 
Ctf{M w ) « -0.12, 

Cgf(M w ) = 0. (6) 

Note that C?nk R in the SM is suppressed compared to C?nk L by m s /m b , because W boson 
couples only to the left-handed fermions. Such terms proportional to m s will be neglected 
in our work by setting m s = whenever they appear. On the other hand, this chirality 
suppression needs not be true in the presence of new physics such as Left-Right symmetric 
(LR) model or in a certain class of supersymmetric models with specific flavor symmetries. 
Such new physics contributions can be parametrized in terms of four complex parameters, 

CT(mw) = Cff (m w ) (£ r -l), 
C8r(m w ) = Ctf{m w ) (6-1), 
C^(m w ) = C^(m w ) 

CT(mw) = Ctf{m w ) £?, (7) 

where £ 7 8 are new complex numbers, whose phases parametrize the effects of the new sources 
of CP violation beyond the KM phase in the SM. The SM case corresponds to £7.8 = 1 and 
£7*8 = 0- R is convenient to define the ratio x as following : 

v = (f - 1W - I) = CfriM w )/Ctf{M w ) 

x-{$8 ij/^ 7 1; as£«(M w )/cffi(M w y K) 

In many interesting cases, this parameter \ is real as assumed in this work. 

Implications of new physics contributions to b —>■ sg have been discussed by various 
group in conjunction with the possible solutions for the discrepancies between theoretical 
expectations and the data on the semileptonic branching ratio of and the missing charms 
in B meson decays, and the unexpectedly large branching ratio for B — > 77 + X s . It has 
been advocated that B(B — > X sg ) ps 10%[~ 50 x B SM (B — > X sg )] can solve these problems 
simultaneously 0. However, this claim is now being challenged by the new measurement 
B(B — > X sg ) < 6.8% (90% CL) 0. In this work, we impose this new experimental data, 
rather than assume that the B(B — > X sg ) is large enough to solve the aformentioned puzzles 
in B decays. 

In the presence of new physics contributions to b — > sj, there should be also generic new 
physics contributions to b —>■ sl + l~ through electromagnetic penguin diagrams. This effect 
will modify the Wilson coefficient Cg of the dim-6 local operator O9 : 
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H eS (b - sll) D # eff (& - s 7 ) - ^ [C 9 9 + C 10 O 10 ] , (9) 

where 

2 2 

° 9 = 4^ ( * L7A) (l7M °' ° 10 = 4^ ( * L7/A) (Ir75/) - (10) 
In the SM, the Wilson coefficients C 9j i 's are given by 

Cl u {M w ) « 2.01, Cf M (Mvr) « 4.55. (11) 

Let us parametrize the new physics contribution to Cg in terms of £ 9 (or x) as following : 

C^{M W ) = C 9 SM (M W )(6 - 1) = Cl M {M w ) x \i 7 - 1), 
, _ C$™(M W )/C$™(M W ) 
X C?£™(M w )/C$™(M w y 1 J 

Since we assume that the new physics modifies only b — > s 7 * and b — > sg, we have C\q(Mw) = 
C^(M^) 0. There is no model-independent relation between £7 and £ 9 , although they are 
generate by the same Feynman diagrams for b — *> s 7 *. In Sec. IV, we will encounter examples 
for both x — and x 7^ in general SUSY models with gluino-mediated flavor changing 
neutral current (FCNC). In principle, there are many more dim-6 local operators that might 
contribute to b — > sl + l~ ||. In the presence of so many new parameters, it is difficult to 
figure out which operators are induced by new physics, since we are afforded only a few 
physical observables, such as B{B — > X a j), B(B — > X s l + l~), Afb(B — > X s l + l~) and the tau 
polarization asymmetry P T in B — > X s t + t~ . Therefore, it would be more meaningful to 
consider the simpler case before we take into account the most general case to figure out 
which operators are significantly affected by new physics. 

Up to now, we considered Or,s(L,R), Og and 0\q relevant to B — > X s 7, X sg , X s l + l~ , assum- 
ing new physics significantly contributes to b — > 57* and b — > sg through dim-5 operators, 
Eqs. (2)-(5). In doing so, five more complex numbers (£7,8(£,i?) > £9) have been introduced. 
If we further assume that the new physics does not induce new operators that are absent 
in the SM, we can drop 0(y i8 )R by setting £ 7)fl = £ 8 fi = 0, thereby reducing the number of 
new parameters characterizing new physics effects into three complex numbers ^7,8,9's (or, 
equivalently £7, x an d X )• Still the number of new parameters are larger than the physical 
observables at our disposal. However, in many interesting cases (and especially SUSY mod- 
els with gluino-mediated b — > s transition that is to be described in Sec. IV), it turns out 
that both x and x are rea l- Therefore, we will assume that both x an d x are rea l hereafter, 
and we are end up with 4 real parameters, which we choose to be |£ 7 |,Im (&),x an d X ■ 
Then we can overconstrain these parameters from the following observables : 

• the branching ratio for B — > X s ~f relative to the SM prediction (i? 7 ) 



2 The Z penguin contribution to b — > sl + l is supressed relative to the photonic penguin by a 
factor of 0(M^/M^), and thus neglected in this work. 
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• the direct CP violation in B — > X s ^j (A b ^, sl = Aq-p) 

• the branching ratio for B — > X sg relative to the SM prediction (R g ) 

• the branching ratio for B — > X s l + l~ relative to the SM prediction (Ru) 

• the forward-backward asymmetry in B — > X s l + l~ (A b f^ sl1 ) 

• the tau polarization asymmetry in B — > X s t + t~ (P t ) 

At this point, it is timely to recall that there have been several works on the model- 
indepedent determination of the Wilson coefficients, CV,(8),9,io from i? 7 and the kinematic 



distributions in Ru [11]]. Our work is different from these previous works in a few 
aspects. First of all, we include the possibility that there is a new physics contribution to Cj 
with a new CP violating phase (Im (£ 7 ) 7^ 0). This necessarily calls for studying the direct 
CP violation in B — > Y s 7 as advocated by Kagan and Neubert || , and invalidates the most 
previous works on the model-independent determination of CV^io's. Secondly we include 
the recent experimental constraint on R g , instead assuming that it can be large enough to 
solve the semileptonic branching ratio problem in B decays. Finally, we assume that the 
new physics does not introduce any new operators with chiralities different from those in 
the SM, and simply modifies the Wilson coefficients of 07,8,9. Thus our analysis does not 
consider the left-right symmetric extension of the SM. 

This paper is organized as follows. In Sec. II, we give basic formulae for the relevant 
physical observables such as _R 7 , Ru, etc. as functions of four real parameters, |£ 7 |, Im (£ 7 ), x 
and x ■ I n Sec. Ill, we present the model-independent numerical analysis for both x = 
and x cases. We show the possible ranges of A CP ,Ru, etc., when we impose the 
experimental data on R 1 and R g . In Sec. IV, we discuss explicit SUSY models with gluino- 
mediated FCNC that enjoy the several assumptions we make in this work. The results of 
this work are summarized in Sec. V. 



II. RELEVANT PHYSICAL OBSERVABLES 

A. B -> V s 7 and B -> X sg 

In the SM, the branching ratios for B — > Xg'j and B — > X sg are obtained including the 
0(a s ) corrections and the nonperturbative effects of b— quark's Fermi motion inside the B 
meson. Relegating the details to the recent works by Kagan and Neubert we show the 
final expressions that will be used in the following : 

^7 = / ( f X y ] \ = 1 + r iW[ Re ^ - 1] +^(x)[|6l 2 - i], (13) 

Rg = n {B (T ^ X y\ = 1 + r 3 (x){Re&) - 1] + r 4 ( X ) [|£ 7 | 2 - 1]. (14) 

-DSMl-D — > y^sg) 

For a real x (and E™ n = 1.95 GeV for the case of i? 7 ), the functions r^s can be approximated 
by§ 
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n(x) ~ 0.46 + 0.020x - 0.0027x 2 , 
r 2 (x) ~ 0.11 + 0.025x + 0.0013x 2 , 
r 3 (x)~0.43x(l-x)+0.50 X , 

r 4 (x) w 0.21 X 2 . (15) 



The recent CLEO data 



£(£ -> X sl ) = (3.15 ± 0.35 stat ± 0.32 syst ± 0.26 model ) x 10" 4 , [12] 

B(B -> X S9 ) < 6.8% (90%C.L.), (16) 

and the SM predictions on these decays ((3.29 ± 0.33) x 1CT 4 ) imply that 

0.77 < Ry < 1.15 (68%C.L.), 

Rg < 6.8/0.2 = 34 (90%C.L.). (17) 
CP violation in the inclusive B — > X s 7 (E™ m = 1.85 Gev) is characterized by CP asymmetry 

a 

A CP = {1.23 lm[C 2 C;\ - 9.52 Im[C 8 C 7 *] + 0.10 Im[C 2 C*]} (%) , 

I w| 

Mx) Mb) (%) j (18) 



A 2 (x)+A 3 ( X ) \&\ 2 + A 4 ( X ) Mb 



where 



A x (x) w0.37-0.18x, 

A 2 (x) w 0.033 - 0.0034x + 0.000085x 2 , 

Mx) w 0.018 + 0.0025x + 0.000085x 2 , 

Mx) w 0.049 + 0.00089x - 0.00017x 2 . (19) 



B. B -> z s z+r 

Now let us consider the decay Z? — > X s l + l~, which occurs through the electroweak 
penguin diagrams and the box diagrams in the SM. If there is a new physics beyond the 
SM, there would be generically dim-6 operators with chiralities different from Og^o shown 
above through the electroweak penguin diagrams and the box diagrams. Morozumi et al. 
considered effects of such new operators (10 operators) on the branching ratio and the 
forward-backward asymmetry (Afb) in 6 sl + l~ ||. In our opinion, it would be more 
meaningful to consider the effects of modified Cj t $ on the decay B — > X s l + l~, since they 
are generically given by dim.-5 local operators. Especially the effects of CV is enhanced by 



3 In this work, the Wilson coefficients without argument represent those at the scale fj, = m^, 
whereas those at the Myy scale are written as Ci(M\y) explicitly. 
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1/s factor in the low s region (see the third line of Eq. ( p0|) below). In any rate, we assume 
that the new physics does not introduce new operators with chiralities different from those 
in the SM, so that we assume that the new physics affects the b — > sl + l~ only through 
modification of the b — > 37*. Therefore, the Wilson coefficients CV^g may change (with a 
new CP- violating phase), and C10 will not be affected at all in our case. 
The differential branching ratio for b — ► sl + l~ is given by || 



dB(b^sl+r) , a 

-1 = B(b^ cez/ 



4vr 2 



1 



f ph (m c /m b )K(m c /m b ) 

y\Cf\ 2 ai(s,m s ,rhi) + |Ci | 2 a 2 (s, rn s , rhi) 

4 \ 
+||C 7 | 2 a 3 (s,m s ,m2) + 12a 4 (s, m 8 , mj) Re {C?Cf }) 



u(s)\l 



Am 2 



(20) 



where all the Wilson coefficients are evaluated at /i = m& by the renormalization group 
equations, s = m^/mf, the function f p h{x) = 1 

space factor for the semileptonic b decays, and the function k(z) defined as 



x 2 + 8a; 6 — x 8 — 24x 4 In x is the phase 



k(z) = 1 



3tt 



( 7r 2_31 )(1 _^ )2 + 3 
v 4 ; 2 



(21) 



is the QCD correction factor thereof. The effective Wilson coefficient C| ff is defined as 

Cf = C 9 f](s)+Y pevt (s) 

= C 9 fj(s) + h{z, s) ( 3 cf + cf + 3Cf + cf + 3Cf + cf ) 

(o) . , M ^ 



1 h(0, s) (cf + 3Cf ) + I (scf + Cf + 3Cf + cf 

2 ^ ' 9 ^ 



(22) 



where s are the Wilson coefficients at // = m;, in the leading logarithmic approximation 



Cf = EM". (j = l,2,...6) 



(23) 



i=l 



with r] = a s (Mw)/ois(fJb), and the numbers a«'s and fc^'s are given in Table XXVII in 
Ref. [I3~|. The functions a's and u are [13 



aci(x,y, z 
a 2 (x,y,z 
a 3 (x,y,z 
a 4 (x,y,z 
uj{s 



2z 



x 



2x 2 + x(l + y) + (1 - yf 

2z 



-2x 2 + x(l + y) + (1 - y) 



+ 



x 



Ax 1 - 5(1 + y)x + (1 - yY 



2z 
1 + — 

a; 

2^ 

1 + — 

x 



-(1 + j/)x 2 - (1 + Uy + y 2 )x + 2(1 +y)(l- y) 2 
(l-y) 2 -(l + y)x 



0s-(l + m s ) 2 ] [s-(l-m s )2], 



(24) 
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with m s = m\jm\ = 0. And fj(s) is given by [[| 
^(s) = 1 + w(s), 

7T 

"00 = - - W " 3 ^ HI -*)- 3 TIT ^y ln(l - s) 
2g(l + a)(l-2a ) lng+ 5 + 93 - 6^ 



3(l-s) 2 (l + 2s) 6(l-s)(l + 2s) 

The function l^> e rt(s) represents the 0(a s ) corrections of the matrix elements, whose explicit 
form can be found at Ref. M. The new physics contributions can induce b — > sqq through 
b — > sg* — > sgg. This will modifies the Wilson coefficients Ci=i-6's, whose effects can be 
seen in the direct CP violation in the B decay amplitude. However these will not affect 
b — > 57 and b — > sl + l~ at the order we are working on. For the realistic prediction, one 
also has to include the long distance contribution through b — > (J/ip^ip ) + s followed by 
(J/ip, ip ) — > . This can be taken into account by adding to the perturbative Y pert (s) the 
resonance contributions [15| : 

with k = —1. To avoid the large contributions from the J/tp and ip resonances, we consider 
the following two regions: the low s region, 1 GeV 2 < s m 2 < 8 GeV 2 for b — > se + e~ case, 
and the high s region, 0.6 < s < 1 for b — > st + t~. 

Using these informations, it is straightforward to evaluate Ru : 

For the decay B — ► X s e + e~, 

B(B -> X s e + e~) = Bi + B 2 \C 9 \ 2 + B 3 Re(C 9 ) + 5 4 Im(C 9 ) 

+ B 5 \C 7 \ 2 + B 6 Re(C;C 9 ) + B 7 Re(C 7 ) + B 8 lm(C 7 ) , (28) 

with 

B x « 1.89, 5 2 « 0.07, 5 3 « 0.19, 5 4 « 0.007, 

fi 5 « 4.07, B 6 « 0.68, 5 7 « 0.87, 5 8 « 0.034 (xl0~ 6 ). (29) 



For the decay B — > X s r + r 



with 



£(£ -> X s r+r-) = D 1 + L> 2 |C 9 | 2 + £> 3 Re(C 9 ) + £> 4 Im(C 9 ) 

+ D 5 \C 7 \ 2 + D 6 Re(C* 7 C 9 ) + D 7 Re(C 7 ) + £> 8 Im(C 7 ), (30) 



Da « 12.6, £) 2 ~ 0.87, L> 3 « -1.51, D 4 « 1.66, 

L> 5 « 6.24, L> 6 » 4.58, L> 7 « -4.21, L> 8 « 4.32 (xl0~ 8 ). (31) 
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Another interesting observable at B factories is the forward-backward asymmetry of the 
lepton in the center of mass frame of the lepton pair : 

__ [Jo 1 rf(cos 6) - $\ d(cos 0)] d 2 B/dsd cos 
FB = [Jo 1 d(cos 0) + d(cos 0)] cPB/dsd cos 



3ogyi - W/j C 10 Re {s[C 9 + y(a)] + 2C 7 } 
{|C 9 + r(s)| 2 ai + C?„a 2 + (4/s)C 7 2 a 3 + 12a JLeC 7 [C g + Y (§)]}' { } 

where 9 is the angle between the positively charged lepton and the B flight direction in 
the rest frame of the dilepton system. For the decay B — > X s e + e~ , the integrated forward- 
backward asymmetry is given by 

, , , Ei + E 2 Re(C 9 ) + E 3 Re(C 7 ) 

Afb ^~ B (B X s e+e~) ' (33) 

where 

E x « 0.21, E 2 ss 0.14, £ 3 « 1.69 (xl0~ 6 ) . (34) 
For the decay B — > X s t + t~ , 

, , s ^ + F 2 Re(C 9 ) + F 3 Re(C 7 ) 

A fb {tt) = — — — — — , (35) 



B(B -> X s r-r- 



where 



Fi « -0.86, F 2 » 1.26, F 3 « 3.66 (xlO^ 8 ). (36) 
The last observable we discuss is the tau polarization asymmetry P T (s) in B — > X s t + t~ 



dB dB 

ds A=-l ds A=+l 

dB i dB 

ds A=-l ~ t ~ ds A=+l 



(37) 



defined as [|16 

A: 

-2u{s)yJ\ - 4m?/sC w Re[(l + 2s) {C 9 + + 6C 7 ] 

~ {\C 9 + F(s)| 2 «! + Cf a 2 + (4/s)C 7 2 a 3 + 12« 4 ReC 7 [C 9 + Y(s)]} 

The integrated tau polarization asymmetry P T can be expressed as 

_ Ti + T 2 Re(C 9 ) + r 3 Re(C 7 ) 

B(B ^ X s r+r-) ' (38) 

where 

Ti » -1.99, T 2 » 2.83, T 3 » 7.31 (xl0~ 8 ) . (39) 

Since B decays into the tau pair probes high m TT (> 3.554 GeV) region, the observable P T 
is sensitive to the deviation of C 9 from their SM values which dominates the B — ► X s t + t~ 
at high s region. 
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III. MODEL-INDEPENDENT ANALYSIS 



Now we are ready to do a model-independent analysis using the formulae obtained in 
the previous section. There are two different cases depending on \ = or not. In principle, 
any new physics contributing to magnetic form factor in b — > 37 may affect the electric form 
factor as well. Therefore one would expect generically x 7^ 0- However this needs not be 
necessarily true as discussed in the next section (the case (i)). So we discuss x — and 
X / fl separately in this section. 

Our strategy is the following : impose the experimental data on i? 7 and R g : 



El : 0.77 < i? 7 < 1.15 as in Ref. p 



• E2 : R g < (6.8%/0.2%) = 34 

For given x an d X 1 these constraints (El) and (E2) determine the allowed region in the 
complex £ 7 plane. Then, in the allowed £7 plane, one can calculate other physical observables, 
Acp, Ru, ApB{b — > sll) and P T . Because the number of observables are greater than the 
number of unknown parameters (one complex number £7 and two real numbers x an d X )> 
one can overconstrain these 4 real parameters. If there is no consistent solution, there would 
be a few possibilities : x and/or x ma y be complex, C w is modified by new physics effects, 
or one has to enlarge the operator basis by including operators with different chiralities from 
those in the SM, as in Ref. [§]. 

Let us first consider the case with x =0. In Fig. [1|, we show the scattered plots of 
various observables as functions of i? 7 for x — 0. The SM case is denoted by a square, 
possible values in our model are represented by dots, whereas the filled circles represent the 
case where there is no new CP violating phase, namely rm£ 7 = 0, but Re^7 7^ 0. Implications 
of these figures are clear. For example, the CP asymmetry in b — > 37 cannot be larger than 
~ ±8% if x — 0) an d Ree can be anywhere between 0.98 to 2.2. For comparison, let us 
discuss the minimal SUGRA model with universal soft mass terms at GUT scale, in which 
typical values of x an d X are X ~ 1 an d X ~ 0.05, x ~ 0.05 respectively JR]. Therefore, 
the predictions in the minimal SUGRA model are very close to the dots in Fig. [l| Namely, 
in the SUGRA case, there are two bands for the possible R ee for a given i? 7 , whereas in our 
case, R ee can be anywhere in between becasue of the presence of a new CP-violating phase 
given by Arctan(Im£ 7 /Re£7). 

In Figs. |2| and|3|, we show similar plots for x — 5 and x = ~~ 5, respectively. This choice of 
X covers a large class of new physics as discussed in Ref. || . Implications of these figures are 
almost the same as Fig. 1, except that there is now rather strong constraint from B — > X sg 
(E2). In this case we can have larger direct CP violation in B — > X s 7 upto 10 — 30%. Also 
the (E2) constraint removes substantial parts of available R ee , Afb(gg) a nd R TT as shown 
in Fig. @ (x = 5), compared to Fig. [I] where the constraint (E2) was not imposed. This 
effect is much more prominent for negative x as shown in Fig. [3] (x = —5). For example, the 
Apsiee) — i? 7 correlation is almost identical to the case with vanishing new phase Im£ 7 = 0. 
From Figs. [1]-^], it is clear that the existence of a new CP violating phase not only can 
generate a large CP asymmetry in b — > 57, but can it also induce quite a lot deviations of 
various observables in b — > sl + l~ for I — e,fi and r. For x — and x = 5, deviations of the 
observables Aq P , R ee , A FB (ee), R TT from their SM values can be large enough that they can 
be clearly observed at future B factories, whereas deviations of other observables Afb{ tt ) 
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and P T from their SM values are rather small that it would be very difficult to measure 
them. For x — ~ 5, only A CP and A FB (ee) shows substantial deviations from the SM values 
because of the (E2) constraint again. If the experimental data on Afb{tt) and P T show 
large deviations from their SM values, it would indicate that x and/or x are complex, or 
some new physics contributes to Cio (with a possibly new CP violating phase), and/or even 
generates Or(s),R and possibly other dimension-6 bsll operators with different chiralities from 
09,io in the SM. 

The nonvanishing x does not affect b — > 37 and b — > sg so that the allowed region in 
the complex £ 7 plane remains the same as before, for a given x- However, it does change 
the observables related with B — ► X s l + l~ , and we show them in Fig. f|for (x,x) — (5, 0.3), 
where we chose x = 0.3 that is typical in the gluino-mediated SUSY models considered in 
the next section. The R ee , R TT , P T dependence on R 1 differ from those in Fig. |], and the 
possible deviations of these observables from their SM values are smaller if x — 0.3. If 
there is no new CP violating phase, the differences are so tiny that one may not be able to 
distinguish two cases in practice. 



The message of this model-independent study is that the previous methods ||- (TT[ has 
to be enlarged to include a new observable Aqp that could be sensitive to a new CP violating 
phase. In the presence of such a new phase, simple correlations among various observables in 
B — > X s j and B — > X s l + l~ (namely, correlations among R 1) Ru, A FB (U) with (/ = e,/x, r) 
and P T that are represented as thick dots in Figs. 1-4 simply disappear, and there is no 
more apparent correlations among these observables. Still one can perform a global analysis 
as before using the formulae given in the previous sections, including the observable Aqp- 
This will provide additional information and one can overconstrain four real parameters, 
Re^7, Im£ 7 , x and x ■ If there is no consistent solution for these four real parameters, one 
has to consider the possibility that x and/or x are a ls° complex. In this case one may be 
able to determine the Wilson coefficients, if one can measure all the observables related with 
B — > X s l + l~ . This task will be possible, only after B factories accumulate the data for the 
first several years. Or one might have to consider the modified 0io and new operators that 
are not possible in the SM. 



IV. SPECIFIC SUSY MODELS WITH GLUINO-MEDIATED FCNC 

A. Models 

In the previous section, we presented model-independent analysis of physics related with 
07,8,9 assuming there is a new CP-violating phase and both x an d X are rea l- m this section, 
we wish to present specific models that satisfy such assumptions. Let us consider the FCNC 
in generalized SUSY models, in which squark mass matrices are nondiagonal in the basis 
where fermion mass matrices are diagonal. In this case, there can be a potentially important 
contributions to the FCNC processes and CP violation that arise from flavor changing f — f — 
g vertices 0. The sources of SUSY FCNC are the nondiagonal (Mf L ) 2 , (M| H ) 2 and (M| R ) 2 . 
Different SUSY breaking models have different patterns/hierarchies for the flavor mixings 
in the squark mass matrices. Since we study the new physics contributions to 07,8,9, the 
Wilson coefficients of the operators already present in the SM in this work, we will consider 
only two cases : (i) the {LR) mixing dominating and (ii) the {LL) mixing dominating cases. 
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There are some models in the literature which fall into these two categories. As discussed 
below, the case (i) does not contribute to C 9 so that £ 9 = 1 (or, x =0). On the other hand, 
the case (ii) contributes to C 9 as well as to £ 7 and £ 8 . Also, there would be generically other 
contributions from H — t,x ~ t an d X° ~ ^fe loops 0. If these loop effects are competent 
with the gluino-mediated loop effects we consider in the following, then our assumption that 
both x and x are rea l would not be true any longer. In the following, we assume that these 
(SUSY) electroweak loops are indeed negligible compared to the gluino mediated FCNC loop 
amplitudes. The latter is enhanced by a s /(Gfm^), as usually assumed. However there is 
a suppression factor in the latter case, the mixing angle in the squark sector given by Tq L 
(or, {S^ll in the mass insertion approximation). Also the heavy squark-gluino loops will 
be suppressed compared to the charged Higgs - top, chargino - stop and neutralino - down 
squarks, unless all the SUSY particles have similar masses so that squark and gluinos are 
not too heavy. So one has to keep in mind that our assumption may break down for too 
small mixing angle in the squark sector or too heavy squark/gluino. With this caveat in 
mind, new physics contributions considered here depend on only one new phase so that x 
and x are rea l; as assumed in the previous section. 

In order to estimate the £7,8,9 in the generalized SUSY models with gluino-mediated 
FCNC, we consider both the vertex mixing (VM) method and the mass insertion approxi- 
mation (MIA). The latter approximation is good, when squarks are almost degenarate. The 
corresponding expressions can be obtained from the former expressions by taking a suitable 
expansion in Am 2 = ((M d ) 2 — m 2 ), where m is a suitable average mass of almost degenerate 
squarks. On the other hand, in the scenario in which the SUSY FCNC and SUSY CP prob- 
lem are solved by decoupling of the (nearly degenerate) first two generation squarks such as 
in the effective SUSY models, there is a large hierarchy between the first two and the third 
squarks so that the MIA is no longer a good approximation. In such case, we have to resort 
to the VM method. 

The full expressions for the Wilson coefficients 6*7.8,9 due to the FCNC gluino exchange 



diagrams are [14 



C* USY (M W ) 



qSUSY 



(M, 



Wj 



X 



°9 



(M w ) 



9V2G F ml\ t 7 tl 



J2 x i( T gl)2I 



vn~ n 



mi) 



J2 X i( F Gl)2I 



(40) 



4 After we submitted this paper, there appeared works which considered these effects in the most 
general MSSM flq|, in the minimal supergravity scenario M] and its modified versions |PQ|. 
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where xj = m~/rn?~ . Tq L and Tq R determine the g — di — dj vertices as follows : 



c = -g s V2(T a ) aP & \(r d GL )j,p L + (r^pj d ja d* 



(41) 



with I,J— 1,2, ...6 and i,j = 1,2,3. They are related with the mixing matrix elements 
diagonalizing the down-squark mass matrix via d~i = {UD)u{d~Li diij'j, UdM^Uq = diagonal, 
with the following identification : 

(IcJ/i = (Uoh, (tfudij = -(Ud)i,j+3. (42) 
The functions /j's are given by JHJ 

x 3 — 6x 2 + 3x + 2 + 6x In x 



fi{x) 
fs(x) 

h{x) 

hix) 



12(x-l) 4 
2x 3 + 3x 2 — Qx + 1 — Qx 2 In x 
12(x - l) 4 : 
x 2 - Ax + 3 + 21nx 



2(x - l) 3 
x 2 — 1 — 2x\nx 



2(x-lf ' 
-llx 3 + 18x 2 - 9x + 2 + 6x 3 lnx 



(43) 



36(x-l) 4 

The corresponding expressions in the MIA is obtained from the above expressions by making 
a Taylor expansion around fh as follows : Xj = (m 2 + Am| ) / m 2 and using the unitarity 

condition for Up. This way one can recover the results in Ref. | I7[| . For completeness, we 
present the resulting expressions below : 



C* USY (M W ) 



S\f2G F m 2 \ t 



(5i) LL M 3 (x) + (5, 

LL \3 



^Mx(x) 



°23 



M 3 (x) +3M 4 (x)) + 



V2G F m 2 X t 
Cl USY {M w ) = (5 d 23 ) P x {x). 

The functions M 1)3 (x) and P\(x) are defined as 

1 + Ax - 5x 2 + (Ax + 2x 2 ) lnx 



L R rrib \3 



M 1 (x)+3M 2 (x) 



(44) 



M x (x) 



Mo(x) = -ar 



2(1 -x) 4 
, 5 - Ax - x 2 + (2 + Ax) In x 
2(1 -xY ' 



M 3 (x) 
M A {x) 
Pi(x) 



-1 + 9x + 9x 2 - 17x 3 + (18x 2 + 6x 3 ) Ins 
12(x - l) 5 : 
— 1 — 9x + 9x 2 + x 3 — 6(x + x 2 ) In x 
6(x - l) 5 ' 
1 - 6x + 18x 2 - 10x 3 - 3x 4 + 12x 3 lnx 
18(x- l) 5 ' 



(45) 
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In order to estimate the £7, x and X > we assume that the (23) mixing is the same order 
of the corresponding CKM matrix element with an unknown new phase <fi ~ 0(1). For 
example, #23 ~ I A t I x e* = AX 2 e** with A = sin# c = 0.22 for both cases (i) and (ii), and 
similarly for Tq L gr . Then it is clear that x and X are rea l i n the MIA both in the cases (i) 
and (ii). In case of the VM approximation, the relevant model is the effective SUSY model 
where only the third family squark can be lighter than ~ 1 TeV so that x±, x-i ~ and we 
may keep only terms proportional to £3 in the summation over I = 1 — 6 in Eqs. (40). Then, 
the x and X are rea l again, as assumed in the previous section. Finally, in the following 
subsection, we will consider only two observables Aqp and R ee for simplicity among several 
observables considered in the previous section. These two observables will be sufficient for 
us to find out the generic features considered in the previous section in the specific SUSY 
models with gluino mediated FCNC. 



B. Case (i) : (LR) insertion 

Let us first discuss the case (i) : (LR) insertion. Since the flavor changing (LR) mixing 
terms are not generated by SUSY breaking in the limit of vanishing Yukawa couplings, 
they are proportional to the corresponding Yukawa couplings. Therefore, the mass insertion 
approximation is always appropriate, and we consider the (LR) insertion only in the MIA. 
From Eqs. (44), one gets 

1 87ra s . . m- Q ^ > 

Cfl\m w ) 9\ / 2G F m 2 X t m b 
= 9 C*f(M w ) |Mx(x) + 3M 2 (x) 
X ~ 8 Ci™(M w ) M,(x) 

X = 0. (46) 

Note that x and x are functions of x only, whereas £7 depends on m, x and also on (^3) lr- 
Therefore, for a fixed x and assuming \8lr\ = \Xt\, one can calculate the Aqp as a function 
of m and (f) with the constraints (El) and (E2). The result is that only x % 1 is consistent 
with the constraints (El) and (E2). As x increases, the contribution to _R 7 and/or R g get(s) 
too large. 

For x = 0.3 and x = 0.8, the allowed range of A C p and R ee as functions of <p are shown 
in Fig. |5| and Fig. ^, respectively, along with the constant fh contours. In the present case 
where the MIA is appropriate, one also has to take into account the constraints on squark 
masses from CDF (m > 230 GeV for x = 1) fH and DO (m > 260 GeV for x = 1) ||| . 



For i^l, one can read off the allowed mass range for the squark mass from the (rrig,rh) 
exclusion plot [p3[| . Roughly speaking, 777 > 300 GeV for x = 0.3 and 777 > 200 GeV for 
x = 3.0. Fig. ^| (a) (x = 0.3 for which x — 1-73, x — 0) indicates that the direct CP 
asymmetry A^~p S7 is in the range < 1.3% for the squark mass m = 400 — 1000 GeV and 
the new CP violating phase = — 0.47T. This aysmmetry is probably too small to be 
observed. But for the same range of 777, and <ft, the R ee can be as large as 2.1 (see Fig. |5| 
(b)). So b —>■ se + e~ is more sensitive to the (LR) mixing than the direct CP asymmetry in 
b — > 57 if x = 0.3. From Fig. |] (a) (x = 0.8 for which x — 5.47 and x = 0); the Aqp is in 
the range 8 — 11% for \cf>\ = 0.2 — 0.357T. It seems that there is a definite lower bound to the 
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Acp, but this is an artifact due to our choice of m < 1 TeV. For heavier m it vanishes very 
slowly (see Fig. |^ (a) and the following paragraph). However if all the squarks (including the 
third family squarks) are heavier than 0(1) TeV, the motivation for the low energy SUSY 
is lost, since the fine tuning problem is reintroduced. Therefore we think that the condition 
fh ^ 1 TeV is a reasonable requirement in the scenarios for the soft SUSY breakings where 
the MIA is valid. With this caveat, the predicted values for Acp are within reach of the 
B factories. The impact on R ee is less striking than the x = 0.3 case, but there is still a 
modest enhancement upto 1.44 of R ee over its SM value which may be also detectable at B 
factories. 

One interesting feature of the (LR) mixing case is that the observables we show in Figs. |5| 
and ^| can probe the effects of very heavy squark masses fh = 400 (710) — 1000 GeV for 
x = 0.3 (0.8). Moreover, the heavier squarks can generate larger Acp, which may be in 
conflict with the naive expectation based on the decoupling of heavy particles in SUSY 
models. However, this is just an artifact of our requirement fh < 1 TeV, as described at the 
end of the previous paragraph. This is because we have fixed x, since the heavier squark 
mass fh for a fixed x implies the heavier gluino mass nig. Therefore the £7 decreases rather 
slowly as fh increases with a fixed x because of the m g factor in the numerator of the second 
term. In Fig. [F] (a) and (b), we plot the direct CP asymmetry Aqp and R ee as functions of 
fh for x = 0.8. We fixed <fi = 0.3tt and 0.57T. If <ft changes its sign, the direct CP asymmetry 
Aqp also changes its sign. We observe that Aqp is maximized around fh = 1 TeV or so. 
The effects of heavy squarks decouple very slowly for Aqp in the {LR) mixing case. On the 
contrary, the effect on R ee is larger for the lighter squark mass as usual. 



C. Case (ii) : (LL) insertion 

Next let us consider the case (ii) : (LL) insertion. In this case, the SUSY breaking terms 
are the main source of the flavor changing {LL) mixing, which are not related with the 
Yukawa couplings in principle. Therefore, the MIA may not be always valid, depending on 
the superparticle spectra. For example, a class of models ||24|| , [ [25]] falls into this case where 
the {LL) mixing dominates. These models [^3] predict that (23) mixing is order of A 2 . 
The mass spectra of the down-squarks in the model p4] are nearly degenerate, whereas in 
the model |25| only the tL,i?,&L, gauginos and the lightest neutral Higgs are relatively light 
compared to ~ 1 TeV. Therefore, one can use the MIA for the first models pl |, whereas 
one has to use the vertex mixing for the second model . 



Below, we will consider the MIA case first. In the mass insertion approximation, 



X 



X 



9 C*f{M w ) \M z {x) + 3M 4 (x) 
8 Cff (M w ) M 3 {x) 
2P,{x) C*f{M w ) 
M 3 {x) Cf M (M, 



23 J LL 



w 



(47) 



In this case we consider two different choices for | (^23)iz, I i n order to compare our results 
with other existing literatures : | (^23) z,i| = |At| and |(^23)ll| = 0(1) f2"7| . As before, 
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one imposes the experimental informations on R 1 and R g , and gets the allowed regions for 
Aqp and R ee for a given phase 0, as well as the direct search limit on the squark mass from 
CDF and DO. For x = 0.3, 1.0 and 3.0, the (x, x) = (7.28, 0.19), (5.25, 0.27) and (3.83, 0.40), 
respectively. Therefore, the overall features of various observables will be close to Fig. £|, 
except that £7 should be fixed to some definite value. 

In case | (^23) ii| = |At|, there are no visible new physics effects on Aqp and R ee , and the 
CP violating dilepton signals can be complimentary to our study. One always have 

|A CP |<1%, |i? ee - 1| < 0.01 (for \(S 23 ) LL \ = |A t |) (48) 

after all the experimental constraints are imposed. However there may be some visible 
deviation in v4pB(ee) as inferred from Fig. H. 



If |(^23)ll| ~ 0(1) as assumed in Ref. [p7|| , then one expects that A C p can be as large as 
±10% to ±15% for ~ ±0.37r for x ~ 1 — 0.3, although the R ee does not change very much 
from its SM value (Figs. § and |9|). If a: gets larger, the Aqp gets smaller and eventually 
becomes undetectable (e.g., Aqp ^ 2% for x = 3, if we impose m > 200 GeV). In Ref. P7] , 
it was noted that this new CP violating phase could result in the CP violation in the decay 
amplitudes for B — > (0, tc°) +K s at the level of (0.1—0.7) of the SM amplitude depending on 
the squark mass. There would be some intrinsic theoretical uncertainties in such estimates 
of nonleptonic exclusive B decays. On the contrary, the direct CP violation in b — > S7 can 
provide independent informations on | (^23) i£ | with less theoretical uncertainties, since we are 
dealing with the inclusive decay rate. In any rate the observable Aqp can play an important 
role in probing a new CP violating phase in B decays if the condition 1(^23) xx| ~ 0(1) is 
met. 

The vertex mixing case can be obtained from Eq. (40) by following identifications : 

(T d \ f (V L ) hJ (for 1 = 1, 2, 3) 
1 I o (otherwise), 

rrM -/ (for 1 = 1, 2, 3) 

1 GR>I ' J ~ \ -(V R )ij (otherwise). 1 yj 

In the (LL) mixing case we consider here, one has (Vr)^ = 5ij. Also, we assume that 
|(Vl)23| — O(0.1) with a new phases of 0(1). This assumption is motivated by a recent 



model by Kaplan et al. [p8 |, which is a SUSY model of flavor based on the single U(l) 
generating the fermion spectra as well as communicating SUSY breaking to the visible sector. 
In this model, only the third generation squarks are lighter than ~ 1 TeV, and the 1st and 
the 2nd generation squarks simply decouple. Therefore, we can keep only the third family 
squarks (/ = 3) in the sum over the squark mass eigenstates, since others are all heavier 
than 0(1) TeV and/or the relations (40), (49) hold. After one imposes the experimental 
informations on _R 7 and R g , one gets the allowed regions for Aqp and R ee for a given phase 
0, as shown in Figs. 10 (a)-(c) for x = 0.3, 1.0 and 3.0. The corresponding values of 
(x, x) are (9.27, 0.16), (7.50, 0.20) and (6.20,0.26), respectively. In Fig. 10, we superposed 
the contours for three different values of 777,3 — 100, 200, 1000 GeV. In the case only the 
third generation squarks are light, the strongest bound on the lighter stop comes from LEP 
experiments ||29|| , and rh 3 = 100 GeV is not excluded yet by LEP experiment. Therefore, 
one expects that Aqp can be as large as 6 % to 12 % for <fi ~ 0.7tt radian for x = 0.3, 1.0 
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and 3.0 respectively, although the R ee does change very little : 1 ^ R ee ^ 1.1. However, 
this large A CP quickly diminishes as m 3 gets heavier, and \A CP \ < 2% for m 3 >~ 200 GeV. 

Therefore, it is very difficult to see the effects of the (LL) insertion in the effective SUSY 
models (for which the VM method is valid) as well as the case of almost degenarate squarks 
(for which MIA is valid) from R ee . In other words, the (LL) insertion can generate a large 
direct CP violation in b — > s'j if there is a new CP violating phase associated with the 
squark mass matrix, (M^ L ) 2 , whereas there can be no significant change in R ee compared 
with the SM case. Also the deviation from the SM diminish very quickly as stop gets heavier. 
Practically speaking, it would be impossible to notice the new physics effects if m > 200 
GeV unless \5ll\ ~ 0(1), for which new physics signal can be visible for heavier stop until 
rh ~ 400 GeV if x is not too large. This is in contrast to the (LR) mixing case, for which the 
new physics effects can increase R ee up to 2.15, and the Acp can be as large as 11 % for fairly 
large stop, rh ~ 1000 GeV and x = 0.8 (see Fig. 6). Although A CP and R ee are not sensitive 
to the (LL) insertion for \5ll\ ~ \, there is another observable which is complementary to 
our observables : namely, CP violating lepton asymmetry in B decays discussed in Ref. |fj6| . 



For larger \5ll\ ~ 0(1), direct CP violations in nonleptonic B decays through AB = 1 
penguin operators can provide additional informations ||27|| . Again, different channels are 



sensitive to different types of new physics, and it will be helpful to study as many modes as 
possible in order to find out new physics signals at B factories. 



V. CONCLUSIONS 

In conclusion, we considered the possible new physics effects on the b —>■ sl + l~ through 
the modified b —>■ 57 vertex. The CP violation in b — > 37 can be very different from the SM 
expectation (Acp(SM) ~ 0), and the branching ratio and A PPj in b — > sl + l~ can be affected 
by the new physics contributing to b — > sj. In particular, the usual model-independent 
extraction of the Wilson coefficients 67,9,10 ma Y be useless in the presence of new physics 
that modifies the CV^ with a new CP- violating phase (namely, Im (£ 7 ) 7^ 0). Therefore, 
not only is the CP asymmetry in b — > 57 a sensitive probe of new physics that might be 
discovered at B— factories, but also it is indispensable for the model-independent analysis 
of b — > sl + l~ . Search for Aqp is clearly warranted at B— factories. 

We also considered specific models which satisfy our assumptions made in the model- 
independent analysis : namely, generalized SUSY models with gluino-mediated FCNC. In 
the case of (LR) mixing, R ee can be enhanced compared to the SM value. Also the direct 
asymmetry Aq P can be as large as 8 — 11% for x = 0.8 and rh = 1 TeV. In this case, the 
direct asymmetry Aqp is sensitive to the heavy squark masses, since the decoupling occurs 
very slowly, beyond rh = 1.1 TeV (see Fig. 7 (a)). Also there is a lower bound on Aqp since 
all the squarks cannot be simultaneously heavier than 0(1) TeV. This is quite an intersting 
feature of the (LR) mixing scenario. In the (LL) mixing case, there is no observable effects 
both for A CP and R ee if = \X t \. But there can be an appreciable amount of A CP upto 
±15%, if \6ll\ — 0(1) in the MIA. In the (LL) mixing with the VM approximation, one 
may be observable Aqp upto 6 — 12% depending on x and the new CP violating phase <fi for 
1(^)23! ~ 0.1 which is the typical values in the model by Kaplan et al. 

Our study is also complimentary to other previous works, e.g., the dilepton asymmetry 
considered by Randall and Su PB| , and the CP asymmetry in the decay amplitudes for 
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nonleptonic B decays considered by Ciuchini et al. ||27|| . It is very important to measure 
various kinds of CP asymmetries at B factories, especially those CP asymmetries which 
(almost) vanish in the SM like the direct asymmetry in b — > 57 and dilepton asymmetry, in 
order to probe new CP violating phase (s) that may be necessary for us to understand the 
baryon number asymmetry of the universe. Different channels may be sensitive to different 
parameter values in new physics, and thus can provide indepedent informations on new 
physics. 



While we were preparing this manuscript, we received a preprint by Chua et al. 



considering the CP-violation in b — > 37 in supersymmetric models. It somewhat overlaps 
with Sec. IV of our present work. But they did not consider the b — > sg constraint, and get 
somewhat larger Acp asymmetry than our work. 
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FIG. 1. Scattered plots for Aqp, R ee , ApB(ee), Rtt, ^fb(tt) and P T as functions of i? 7 for 
X = X = 0- The SM values are marked as a filled square. 
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FIG. 2. The same as Fig. 1 for x = 5 and % = 0. The SM values are marked as a filled square. 
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FIG. 3. The same as Fig. 1 for x = —5 and x = 0. The SM values are marked as a filled square. 
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FIG. 4. The same as Fig. 1 for x = 5 and x = 0.3. The SM values are marked as a filled square. 



24 





25 




26 





500 1000 1500 2000 2500 3000 3500 

(b) ™ ( GeV ) 



FIG. 7. (a) A°^ sy and (b) R ee as functions of the squark mass rh in the (LR) mixing case with 
x = 0.8. 



27 




-0.5 



a 



2 2.5 

^(tt rad) 



02 



14 
1.12 

1.1 
1 .08 
1 .06 
1 .04 
1 .02 
1 

0.98 



-0.5 



250 



400 



nun 1000 







0.5 



.5 



2 2.5 

^(tt rad) 



FIG. 8. Possible values of Acp and R ee as functions of the new phase 4> in the (LL) mixing 
case with x = 0.3 (i.e., x = 7.28 and % = 0.19). We assume that | (^23)^^! = 1- 
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FIG. 9. Possible values of Acp and R ee as functions of the new phase 4> in the (LL) mixing 
case with x = 1 (i.e., x = 5.25 and %' = 0.27). We assume that | (#23)1,1, | = 1- 



29 



o 

< 



5 _- 
~- 

-5 - 

Z L 

-0.5 



1 00 
200 
000 



^,rT1,T,X,.-..i,,. ., 








0.5 



1 

a 



.5 



2 2.5 

^(tt rad) 



CJ 

< 



-10 - 



-0.5 





1 00 

200 

1 000 


- I I I I I I I I I I I 




J.I l U I mi | (m i lii i Hi" u: ■ uj ■ uj 

i i i i i i i 







0.5 



b 



.5 



2 2.5 

^(tt rad) 



^ 10 - 

cj Z 

< - 

-10 - 



-0.5 







0.5 



1 00 
200 
1 000 




.5 



2 2.5 

^(tt rad) 



FIG. 10. The Aqp contours in the (to, 4>) plane for (a) x = 0.3, (b) x = 1 and (c) x = 3 in the 
(LL) insertion case using the vertex mixing method. 
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